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Abstract. In this paper, we propose an approach to controller synthesis for a class of constrained 
nonlinear systems. It is based on the use of a hybridization, that is a hybrid abstraction of the 
nonlinear dynamics. This abstraction is defined on a triangulation of the state-space where on 
each simplex of the triangulation, the nonlinear dynamics is conservatively approximated by an 
affine system subject to disturbances. Except for the disturbances, this hybridization can be seen 
as a piecewise affine hybrid system on simplices for which appealing control synthesis techniques 
have been developed in the past decade. We extend these techniques to handle systems subject to 
disturbances by synthesizing and coordinating local robust affine controllers defined on the simplices 
of the triangulation. We show that the resulting hybrid controller can be used to control successfully 
the original constrained nonlinear system. Our approach, though conservative, can be fully automated 
and is computationally tractable. To show its effectiveness in practical applications, we apply our 
method to control a pendulum mounted on a cart. 



1. Introduction 

Abstraction based approaches to control synthesis have become quite popular (see [20] and the refer- 
ences therein). These consist in computing a simple conservative approximation (i.e. the abstraction) 
of the system dynamics. A controller synthesized using the abstraction can then be used to control 
the original system. Most of these approaches use discrete abstractions since these allow to leverage 
controller synthesis techniques developed in the area of discrete event systems. The use of hybrid 
abstractions, called hybridizations, though common for verification purpose [131 El US E] , has been 
much less explored in the context of controller synthesis. 

In this paper, we present an approach to control synthesis for a class of constrained nonlinear systems 
based on the use of a hybridization. The abstraction is defined on a partition of the state-space in 
simplices (i.e. a triangulation). In each of these simplices, the nonlinear dynamics is approximated 
by an affine system with disturbances where the disturbances account for the approximation errors. 
The resulting abstraction can be seen as a piecewise affine hybrid system on simplices subject to 
disturbances. For such systems, without disturbances, automated control synthesis techniques have 
been developed by Habets and van Schuppen in a series of papers [HI [121 EH] (see also the related 
work in [H [TH EH]). We extend these techniques to handle the systems subject to disturbances by 
synthesizing local robust affine controllers defined on the simplices of the triangulation. We show 
that by coordinating these local controllers it is then possible to control successfully the original 
constrained nonlinear system. A preliminary version of these results was presented in [9]. The 
technical content in the present version has been improved: the class of considered abstractions is 
more general (the disturbance is constrained in an arbitrary compact convex set); we synthesize 
hybrid controllers instead of piecewise smooth controllers, this allows us to prove the correctness 
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of the synthesized controller more rigorously; finally, a more challenging application example is 
considered. 

The paper is organized as follows. In the next section, we formulate the control problem under 
consideration. In section 3, we introduce the notion of hybridizations and show how these can be 
computed for nonlinear systems. In section 4, we extend the techniques of |X0^ [19] for computing 
local robust affine controllers on simplices. In section 5, we show how to use these local controllers 
and the hybridization in order to synthesize a controller for the original system. Finally, in section 6, 
we show the effectiveness of our approach by applying it to control a pendulum mounted on a cart. 



2. Problem Formulation 



In this paper, we consider a constrained nonlinear system of the form: 

*(i) = /(x(t)) + 0(x(t))u(t), 



x(i) G X, u(t) G U 

where the state domain X C M. n is the union of a finite number of compact polytopes and the 
control domain U C M m is a compact polytope. We shall assume that the maps / : X — > M. n and 
g : X — > W ixrn are of class C 2 and C 1 respectively. Let Xj C X and Xt Q X be compact polytopes, 
specifying a set of initial states and a set of target states, respectively. 

Definition 2.1. A hybrid controller for system S is a tuple C = (Q, E, D, Q, %), where Q is a finite 
set of modes; E C Q x Q is the set of edges; V = {D q C X\ q G Q} is the set of domains with D q 
a compact set for all q G Q; Q = {G e C D q n -D q '| e = (q, (/) G £"} is the set of guards with G e a 
compact set for all e G E; % = : -Dq — >• J7| q G Q} is the set of continuous controllers with h q a 
Lipschitz continuous map for all q G Q. 

The system £ controlled by C is a hybrid system, denoted S||C whose evolution is described by 
executions [15]. An execution of S||C is a (finite or infinite) sequence of triples (r^, qi,^-i)^L consisting 
of time intervals Tj = (if iV is finite we can have r^r = [£tv,+oo)), modes qi G Q, and 

differentiable maps Xj : — )■ A such that to = 0, xo(to) G Xj n -D go ; for all i = 0, . . . , N with 
i < +oo, for all t G 

xi(t) G Z\. and Xi (t) = /(x;(t)) + g(xi(t))h qi (xi(t)); 

and for all i < A, 

(9t,g»+i) G -E, and Xi(t i+ i) = x i+1 (ii + i) G G fe>9i+l) . 
The execution is said infinite if A is infinite or r^r = [tjv,+oo), otherwise it is said finite; it is said 
Zeno if A is infinite and limj_-. +00 ij is finite. S| \C is said to be non-blocking if all finite executions can 
be extended into infinite ones; it is said non-Zeno if it does not have Zeno executions. We consider 
the following synthesis problem: 

Problem 2.2. Synthesize a hybrid controller C, such that Xj C [j q& nE) q and the hybrid system 
S||C is non-blocking, non-Zeno and all its infinite executions (Tj,gj,Xj)^ satisfy 
3i* G N, < i* < A, Vi = i*, . . . , A, Vi G Xj(t) G A T . 



Let us remark that since for all g G Q, -D„ G A, all infinite executions of S||C also satisfy for all 



i = 0, . . . , A, for all t G Tj, Xj(t) G A. Our approach to Problem 2.2 consists of two main ingredients: 
a hybridization of the nonlinear system E and robust continuous controllers defined on simplices. 
These are presented in the following sections. 
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Notations. In the following, for two vectors vi, V2, v\ -V2 denotes their scalar product. For a vector v, 
\\v\\ = y/v ■ v is its Euclidean norm. For a matrix A, \\A\\ denotes its norm induced by the Euclidean 
norm for vectors, and p{A) denotes its spectral radius. 

3. Hybridization 

A hybridization is a hybrid abstraction of a continuous dynamical system pQ. Conservativeness 
of the approximation is ensured by the introduction of disturbances. In this paper, we shall use 
hybridizations for control synthesis. We first describe the hybridization principle. 

Definition 3.1. S = {S p \ p G P} is a triangulation of the domain X if the following conditions 
hold: for all S p G S, S p is a full dimensional simplex of W 1 (i.e. S p is the convex hull of n + 1 affinely 
independent points); for all S p , S p i G S, their intersection is the convex hull of their common vertices 
or empty; \J peP S p = X. 

We further assume that S contains a triangulation of Xt, St Q S. The main idea of the hybridization 
principle consists in approximating the dynamics of E locally, in each simplex S p 6 5, by an affine 
dynamics with disturbances of the form: 

f x(t) = A p x(t) + B p u(t) + a p + w(t), 
\ u(t) G U, w(t) G W p 

where the disturbance set W p is assumed to be a non empty convex compact set; let us define 
S' = {(S p , A p , Bp, a p , W p )\ p G P}. With the exception of the disturbance, the system £' can be seen 
as a piecewise affine hybrid systems on simplices for which control synthesis techniques have been 
developed in [10]. The disturbance is added to compensate the approximation error and guarantees 
the conservativeness of the approximation: 

Definition 3.2. S' is a hybridization of system E if and only if for all p G P, 
(3.1) Vx G S p , Vu G U, (/(x) + g(x)u — A p x — B p u — a p ) G W p . 

We now propose a method for the computation of a hybridization. We do not discuss the computation 
of a triangulation S of the domain X. This is a well studied problem in computational geometry 
for which efficient algorithms exist, at least in low dimensional spaces (see e.g. p2]). I n higher 
dimensional spaces, provided the domain can be partionned in hypercubes, a triangulation of X can 
be obtained using a simple triangulation of each hypercube as shown in [lj. It is to be noted that 
the size of the simplices of the triangulation can generally be made arbitrary small. 

We now focus on the computation of the affine dynamics in simplex S p G S. Let {f P) o, • • • , v Ptn } 
denote the n + 1 vertices of the simplex. Essentially, the map f(x) is approximated by the affine 
vector field A p x + a p while the map g(x) is approximated by the constant matrix B p . A p and a p can 
be obtained by interpolation of / at the vertices of the simplex: 

Mi G {0, . . . , n}, ApVp^i + a p = f(v p>i ). 

Since v Pt o, . . . , v PtTl are affinely independent, this condition uniquely determines A p and a p which 
can be computed by solving a set of linear equations. Then, the matrix B p can be chosen as the 
barycentric combination of the value of g at the vertices of the simplex: 

1 n 

i=0 
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It remains to compute the disturbance set W p such that equation (3.1) holds. For / £ {1, . . . , n}, let 
fi(x) denote the l-th component of the vector /(x), gi,{x) denote the Z-th line I of the matrix g(x), 
a p j denote the l-th element of a p , A p j, and B p i, denote the l-th lines of A p and B p . The set of 
disturbances W p can be chosen as W p = [/jl ^ Mp,i] x ' ' ' x \p n > ~fi P ,n] where for I £ {1, . . . , n}, 

a = min (f t (x) + gi,(x)u - A pl ,x - B pl ,u - a p r) 



and 



H i = max (fi(x) + gi.(x)u - A p i,x - B p i,u - a p i) 

x£S p ,u£U 



Then, it is clear that equation (3.1) holds. The following result shows that for the hybridization 
defined above, the disturbance set W p can be made arbitrarily small provided the triangulation of 
the domain is sufficiently fine. 

Proposition 3.3. Let £' = {(S p , A p , B p , a p , W p )\ p £ P} be the hybridization defined above. Let 
H(fi) denote the Hessian matrix of fi and J{gi,) denote the Jacobian matrix of gi,. Let 5 P denote 
the diameter of the simplex S p : S p = max{||a; — x'\\\ x,x' £ S p }. Then, W p C [— fj, Pl i, (J,p,l\ x • • • x 
[-Mp,n, fJy,n] where for I £ {1, . . . , n}, 

H P d = a P ,l8p + /3 p ,iS p me^ \\u\\ 
with a p j = -maxp(f/(/i)(x)) and /3 p j = max\\J(gi,)(x)\\. 

Proof. It is clearly sufficient to show that for all x £ S p , for all I £ {1, . . . , n} 

\fi(x) - A P)i ,x - a p j\ < a p j5 p and \\gi.(x) - B p>l ,\\ < P p ,i6 p . 



A proof of the first inequality can be found in [21|. The second inequality is derived as follows. From 
the mean value theorem for vector- valued function we have for all x±, x<± in S p 

\\9i»(xi) -9im(x 2 )\\ < max||J(^.)(x)||||xi - x 2 \\ < m&x.\\J(gi,)(x)\\6 p . 

x£jp x(Eijp 

Then, it follows that for all x £ S p , 

Y n 1 n 

\\gi.(x) - B p i.\\ < \\gi.(x) ^5^^.(^,011 < — r~r Y) \\9i»(x) - 9i»(v p ,i)\\ 

n + 1 n + 1 z — ' 

i=0 i=0 

1 n 

- rxrS^l^^^l^p = ™II J (^»)( X )IIV 

□ 

We can see that the bound on the disturbance depends linearly in the diameter of the simplex. This 
means that to reduce the bounds by factor 1/2 it is necessary to consider 2 n times more simplices 
in the trinagulation. Hopefully, as we shall see on a practical example, it is not always necessary to 
consider very fine partitions for our approach to be successful. 



The rest of the paper is devoted to solving Problem |2.2| Our approach is inspired by the work 
presented in [10J for piecewise affine hybrid systems on simplices. 
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We need to extend several techniques developed for the class of afflne systems on simplices to the 
class of systems with disturbances. The results are stated without proofs which are straightforward 
adaptations of the proofs in [TOj, [19] and can be found in [9] . Let S be a simplex of W 1 , we denote 
vq, . . . ,v n and Fq, . . . , F n the vertices and the facets of S with the convention that Fj is the facet 
opposite to vertex Vj. mo, . . . , m n denote the outward unit normal vectors of the facets of S. 

4.1. AfRne systems with disturbances. We consider the following autonomous affine system 
with disturbances : 

(4.1) x(i) = Ax(t) + a + w(t), x(i) G R n , w(i) G W 

where the disturbance set W is a non empty convex compact set. It is sufficient to assume that the 



disturbance w is a continuous map. We say that a trajectory x of (4.1), starting in S, exits S at 
time T > 0, if there exists e > such that 

\ft G [0,7], x(i) G S and Vt G (T, T + e), x(i) ^ 5. 

We shall say that a facet -Fj is blocked if 

Vx G Fi, rn,- • (Ax + a) < — maxm,- • <i. 
J J d&D J 

The fact that a facet is blocked can be characterized using only the value of the vector field at the 
vertices: 

Proposition 4.1. The facet Fj is blocked if and only if 

Vi G {0, . . . , n}, i ^ j, rrij ■ (Avi + a) < — max rrij ■ d. 

The following proposition characterizes the points where the trajectories may exit the simplex. 



Proposition 4.2. // a trajectory x of (4-1), starting in S, exits S at time T, then x(T) belongs to 
a facet of S that is not blocked. 



We now give necessary and sufficient conditions such that all the trajectories of (4.1), starting in S, 
exit S in finite time. 



Proposition 4.3. All the trajectories of (4-D starting in S, exit S in finite time if and only if 



3£ G R n , Vx G S, £ ■ (Ax + a) > - min £ • w. 

The previous property can again be characterized using only the value of the affine vector field at 
the vertices of the simplex: 



Proposition 4.4. All the trajectories of (4-D starting in S, exit S in finite time if and only if 



3£ G M. n , Vi G {0, . . . , n}, £ • (A Vi + a) > - min £ • w. 
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4.2. Robust controller synthesis. Let us now consider an affine control systems with disturbances 
on simplices of the form: 

x(i) = Ax(i) + Bu(t) + a + w(i), x(t) G R n , u(i) G f7, w(t) G W 

where the control domain U C M p is a convex compact polytope and the disturbance set W is a 
non empty convex compact set. The disturbance w is assumed to be continuous. We consider the 
following control problem: 

Problem 4.5. Consider a subset of indices J C {0, ...,n}, and the associated subset of facets 
T = {Fj\, j G iJ}, design an affine feedback controller h : W 1 — > M m , h(x) = Kx + k, such that for 
all x G S, h(x) G U and for the autonomous affine system with disturbances 

(4.2) x(t) = Ax(t) + Bh(x(t)) + a + w(t), x(i) G M n , w(i) G 



all facets that are not in T are blocked and all trajectories of (4.2) starting in S exits S in finite 
time. 

We denote Uq, • • • , u n G U the values of the controller at the vertices of S: 

Ui = h(vi) = Kvi + k, i G {0, . . . , n}. 

Since Vq, . . . , v n are affinely independent, Uq, ■ ■ ■ ,u n uniquely determine the matrix K and the vector 
k. Moreover since no, . . . , u n G U and for all x G S, h{x) is a convex combination of uq, . . . ,u n it 
follows that h(x) G U, for all x G S. At the vertices of the simplex S, the value of the vector field of 



( 4.2 ) is given by Avi + Bh{vi) + a = Avi + Bui + a. In the following we characterize suitable values of 



UQ,...,u n with the understanding that these allows the computation of the affine controller h. The 



following result is a direct consequence of Propositions 4.1 and |4.4| 

Proposition 4.6. For i G {0, . . . , n}, let us consider the following convex polytopes: 



rrij ■ (Avi + Bu + a) < — maxmj • d, Vj G {0, . . . , n} \ {jj U {i}) 



An affine feedback controller h solves Problem 4-5 if and only if for all i G {0, . . . , n}, m G Ui and 



e 1" , Vi G {0, . . . , n}, (Avi + Bm + a) > - min £ • w. 

Let us discuss the computation of input values uo, . . . ,u n satisfying the conditions of the previous 
proposition. We denote by Ui the set of vertices of Ui. 



Proposition 4.7. There exists an affine controller h solving Problem 4-5 if and only if there exists 
uo G Uo, . . . , u n G U n such that 

3£ Vi G {0, . . . , n}, £ • (Avi + Bu { + a) > - min £ • w. 



A controller solving Problem 4.5 can therefore be synthesized by computing the vertices of the 



polytopes Uq, . . . , U n and then looking for a suitable combination of these vertices. 



5. Nonlinear Control Synthesis 



We now describe a solution to Problem |2.2| based on the results described in the previous sections. 
Let £' be a hybridization of E. Our approach is based on the coordination of robust affine controllers 
defined on the simplices of the triangulation S. 
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5.1. Invariance controller. We start by synthesizing a controller in order to render the set Xt 
of target states invariant. For all simplices S p G St of the triangulation of Xt, we denote by 
vertices(Sp) its set of vertices and by external- f acets(S p , Xt) the set of facets of S p that are included 
in the boundary of the polytope Xt- For a given facet F G external-facets(S p , Xt), we denote by 
rriF its outward unit normal vector and by vertices(F) its set of vertices. 

Problem 5.1 (Invariance control). Synthesize a Lipschitz continuous controller Kt '■ Xt —> U, such 
that for all xq G Xt, the solution of 

(5.1) x(i) = /(x(t)) + 0(x(t))Mx(t)), x (0) = ^0 

is defined for all i G M + . 

Let us remark that since Ht is defined only on Xt, this implies that x(i) G for all t G M + . We 
will search the controller hr as a continuous piecewise affine map defined on the triangulation St 
(i.e. Ht is affine on each simplex of the triangulation). Ht is uniquely determined by its value at the 
vertices of the simplices of St- If all these values belong to U then so does hx(x), for all x G Xt- 
Also, it can be shown [8] that Ht, defined this way, is Lispchitz on Xt- Then, from the standard 
characterization of invariant sets [2] we have the following result: 



Lemma 5.2. Ht '■ Xt — > U solves Problem 5.1 if and only if for all S p G St, for all F G 



external-facets(S p , Xt) , for all x G F, txif ■ (f(x) + g(x)hT(x)) < 0. 

We now characterize the values of Ht at the vertices of the simplices of St that result in controllers 



solving Problem 5.1 



Proposition 5.3. If for all S p G St, for all F G external-facets(S p , Xt) , for all v G vertices(F) 

(5.2) rap ■ (A p v + B p hT(v) + a p ) < — max tuf ■ w 

uieWp 



then hT solves Problem 5.1 



Proof. Let S p G St, F G external- f acets(S p , Xt), and x G F, then 

niF ■ (f(x) + g(x)h T (x)) = m F ■ (A p x + B p h T (x) + a p ) 

+m F ■ (f(x) + g(x)h T (x) - A v x - B p h T (x) - a p ). 

Since £' is a hybridization of X, and since h(x) G U, it follows from Definition |3.2| that (f(x) + 
g(x)hT(x) — A p x — B p hT(x) — a p ) G W p . Therefore, 

m F ■ (f{x) + 9(x)hT(x)) < rap ■ {A p x + B p hT(x) + a p ) + max rap ■ w. 



Since A p x + B p hT(x) + a p is an affine map on F, we have by convexity and by (5.2) that 

mp ■ (A p x + B p hp(x) + a p ) + max mp ■ w < 0. 



Then, from Lemma 5.2, we have that the controller hp solves Problem 5.1 □ 



Let us remark that the computation of the controller hT involves only finding values of satisfying 
a set of linear inequalities, at the vertices of the triangulation. This can be done efficiently using 
linear programming. For simplicity of the presentation, we assumed that we render the whole target 



set Xt invariant; actually, as far as Problem 2.2 is concerned, it is sufficient to render a subset 



Xlp C Xt invariant. This can be done a similar way. 
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5.2. Reachability controller. We now describe how to synthesize a hybrid controller solving Prob- 



lem 2.2 The proposed controller essentially drives the trajectories of the system through a sequence 
of simplices ending in Xt- This is done by coordinating robust affine controllers defined in section [4~2{ 



Let S' be a subset of the triangulation S and let S p G S, we denote by common- f acets(S p , S') the 
subset of facets of the simplex S p that are also facets of a simplex in S' . We denote by adjacent{S') the 
subset of simplices that are adjacent to a simplex in S': S p G adjacent{S') if and only if S p G S\S' 
and common- facets(S p ,S') ^ 0. The synthesis of the hybrid controller C can be done using the 
following algorithm: 



Algorithm 5.4. Let Q :- 

Let k := 0, 5_i := 0, S :■- 
While S k / 



{<?t}, D qT 
St- 



Xt, h qT solve Problem 5.1 



E :-- 



S 



fc+i 



For all S p 6 adjacent(Sk 



if there exists an affine feedback controller /i p which solves 
(S p , common- facets(S p ,Sk), A Pl Bp,a p ,W p ), then let 



S^, /ip solution of Problem 4.5 and: 



Problem [45] with (5, T, A, B, a, W) = 
cSfc+i := S k+ i U {S p }, Q := Q U {p}, £> p 

— For S'p' G Sk\ST, if common- facets(S p ,{S p >}) / then E 1 := i?U and G( PjP /\ := 
common-facets(S p , {S p i}), 

— If common-facets(S p , St) / then E 1 := i?U{(p, <?t)}, C(p,ijt) := common-facets(S p , St) ■ 
• k := k + 1. 

Return C = (Q,E,V,g,H) where £> = {£> g | gGQ},e = {G e | ee£},?/ = g G Q}. 



The algorithm essentially applies dynamic programming: Sk consists of the simplices from which we 
can drive the system to Xt, using robust affine controllers, by following a sequence of at most k 
simplices; at each iteration we add to Sk+i the simplices in which we can drive all the trajectories to 



facets that are shared with some simplices in Sk- Let us remark that the Algorithm 5.4 necessarily 
stops after K < |«S| iterations and each iteration involves solving a number of linear programs that 
is at most linear in the number of simplices. Then, the worst-case time complexity of the algorithm 
is quadratic in the number of simplices. However, it is to be noted that the number of simplices 
in the triangulation is exponential in the dimension of the state-space. The algorithm is therefore 
exponential in the dimension of the system. This is expected for an approach based on abstractions 
defined on a partition of the state-space. 



Theorem 5.5. Let C be computed by Algorithm 5.4' if Xj C {J q ^qDq then C solves Problem 2.2 



Proof. Let us consider a finite execution (T^g^x, 



)£o of n\c. 



If q N 



then since h qT solves 



Problem 5.1 x^v can be extended on [tjsr, +oo). If / qT, then gjy G P, let k < K be the smallest 
index such that S„ N G Sk- Since X/ is a hybridization of S, the continuous dynamics in mode q^ is: 



x(t) = /(x(t))+ 5 (x(t))V(x(t)) 

= A qN x(t) + B qN h qN {x(t)) + a qN + w(t), w(i) G W qN 

for an obvious particular value of the disturbance w(t). By construction of h qN , all trajectories of the 
linear system above exits S qN in finite time and by Proposition 4.2 , xat can be extended on [t^, t' N+1 ] 
where xjy(t' N+1 ) belongs to a facet F qN of S qN that is not blocked for the linear system above. By 
construction of the guards, there exists an edge e = (qN, qN+i), such that F qN C G e . Hence, the finite 
execution can be extended into an execution of the form (rj, Xj)^ 1 . Moreover, either qN+i = qr 
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or we have S qN+1 £ Sk-i- From the previous discussion, and using a simple induction, we can show 
that every finite executions of T,\\C can be extended into infinite ones, by taking a finite number of 
discrete transitions (at most K) before reaching the final mode qr in which continuous evolution 
is enabled until time goes to infinity. Hence, S||C is non-blocking. Since all executions have only 
a finite number of discrete transitions, is necessarily non-Zeno. Finally, all executions of S||C 



reach the final mode qx and stays there forever, since D„ T = Xt, Problem 2.2 is solved. □ 



Let us remark that our approach for solving Problem 2.2 is clearly conservative. Our algorithm 
may fail to solve the problem even though a suitable controller exists. There are several sources of 
conservatism. The first one is due to the use of a hybridization. Proposition |3 . 3| suggests that this can 
be reduced by using a finer triangulation at the price of an increased computational effort. The other 
sources of conservatism are inherent to the approach developed in [10]. However, this conservatism 
allows us to synthesize controllers that are correct by design by a fully automated method which is 
computationally effective as shown in the following section. 

6. Example 

In this section, we apply our approach for controlling a pendulum mounted on a cart. This system has 
been considered in [18] for illustrating controller synthesis using discrete abstractions. The dynamics 
of the system is described by 

xi (t) = x 2 (t) 

X2(i) = — sin(xi(t)) — u(t) cos(xi(i)) , u(i) € [— a, a] 

where xi(i) is the angle between the pendulum and the downward vertical, x 2 (i) is the angular 
velocity. The control input u(t) is the acceleration of the cart and is assumed to be bounded in 
absolute value by a. The goal is to bring the pendulum from downward vertical position (xi(i) = 0) 



to upward vertical position (xi(t) = %) and to keep it there. This can be formulated as Problem 2.2 
where X = [0, 27r] x [— it, ft], Xj = {(0,0)} and Xt is a neighborhood of (vr,0). Let us remark that 
since xi(t) represents an angle then if the trajectory exits the domain X with xi(i) = 2ir, it re-enters 
immediately with x\(t) = and the same velocity. This is important when defining the adjacency 
relations between simplices of the triangulation of X. 

We use a regular triangulation defined on a partition of X into squares of length tt/N where iV G N 
is a parameter. Each of these squares is then partionned into two triangles: this can be done using 
a line of slope 1 or —1 giving two possible orientations for the simplices. The target set Xt is given 
by 

The orientation of the simplices is left free at the beginning and is determined on the fly by Al- 



gorithm 5.4 The hybridization is also computed on the fly. Let us remark that only the sec- 
ond component of the vector field is nonlinear, therefore the disturbance set W p is of the form 
W p = {0} x 2 , where the bounds given in Section |3j can be computed explicitly. When com- 
puting the continuous controllers, to choose among the possible values at the vertices of the simplex, 
we try to minimize the variation of the control map within the simplex. In particular, if a constant 
controller solves the problem then such a controller will be used. 



Algorithm |5.4| was implemented in Matlab, and ran on a standard desktop equipped with CPU 
Pentium 4 (3.20 GHz). We report in Table [l] performances of the algorithm for varying number of 
simplices and set of inputs. We can check that the computational costs increases as the constraints on 
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u ^^^^ 


10, 800 


20, 3200 


30, 7200 


40, 12800 


[-2,2] 


6 - No 


176 - Yes 


591 - Yes 


1263 - Yes 


[-3,3] 


6 - No 


102 - Yes 


330 - Yes 


631 - Yes 


[-4,4] 


14 - Yes 


90 - Yes 


251 - Yes 


531 - Yes 



Table 1 . Run times in seconds and success of the synthesis for varying number of 
simplices and set of inputs. 




the inputs become tighter. Also, the complexity with respect to the number of simplices, estimated 
experimentally on successful cases is roughly in 0(|5| 3//2 ) which is polynomial but better than the 
theoretical worst-case complexity which is quadratic. Compared to the approach presented in |18j . 
where the same example is considered, the approach seems to have quite similar performances: 
reported computation times are of the same order and the estimated complexity appears to be 
polynomial in the number of discrete states as well. 

In Figure [TJ we show the controller obtained for U = [—4,4] and N = 10 (i.e. |<S| = 800), as well 
as the trajectory of the controlled system. The pendulum starts balancing in the anti-clockwise 
direction to accumulate some energy, then it moves in the clockwise direction until it reaches the 
target. This nontrivial example shows that our approach can be used for automated controller 
synthesis for nonlinear systems. 



SYNTHESIS FOR CONSTRAINED NONLINEAR SYSTEMS 
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7. Conclusion 

In this paper, we presented an algorithmic approach to controller synthesis for constrained nonlinear 
systems. Our technique is based on two main ingredients, namely a hybridization and robust affine 
controllers on simplices. Though conservative, our method can be fully automated and we showed 
that it is effective on an example. Our method should probably be reserved to small-dimensional 
systems as the number of simplices in the triangulation explodes when the dimension grows. 

There are several possible extensions for this work. First, instead of using piecewise affine hybridiza- 
tion defined on a triangulation, one could use piecewise multi-affine hybridization [1] defined on a 
partition of the state domain in hypercubes as an extension of |3j would allow us to synthesize local 
controllers on hypercubes. Second, by extending the approach presented in [TJ [T3], one could solve 
more complex control problems such as those specified in linear temporal logic (LTL). Finally, the 
class of nonlinear system we consider could be extended, by considering systems of the form 

x(t) = /(x(t),u(t)), x(i) G X, u(t) G U. 

Following (6j, a piecewise affine hybridization could be computed by triangulating the domain X x U. 
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